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For the evolution of density fluctuation in nonlinear cosmological dynamics, adhesion approxi- 
mation (AA) is proposed as a phenomenological model, which is especially useful for describing 
nonlinear evolution. However, the origin of the artificial viscosity in AA is not clarified. Recently, 
Buchert and Dommguez report if the velocity dispersion of the dust fluid is regarded as isotropic, 
it works on a principle similar to viscosity or effective pressure, and they consider isotropic velocity 
dispersion as the origin of the artificial viscosity in AA. They name their model the Euler-Jeans- 
Newton (EJN) model. In this paper, we focus on the velocity distribution in AA and the EJN 
model and examine the time evolution in both models. We find the behavior of AA differs from 
that of the EJN model, i.e., although the peculiar velocity in the EJN model oscillates, that in 
AA is monotonically decelerated due to viscosity without oscillation. Therefore it is hard to regard 
viscosity in AA as effective pressure in the EJN model. 
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I. INTRODUCTION 

The Lagrangian description for the cosmological fluid can be usefully applied to the structure formation scenario. 
O^l ' This description provides a relatively accurate model even in a quasi-linear regime. Zel'dovich pj proposed a linear 
Lagrangian approximation for dust fluid. This approximation is called the Zel'dovich approximation (ZA) p], 0, 
El B . B ijB. El [M 0, E| • ZA describes the evolution of density fluctuation better than the Eulerian approximation 
fllL fia . fl3L fl^ | . Although ZA gives an accurate description until the quasi-linear stage, ZA cannot describe the model 
after the formation of caustics. In ZA, even after the formation of caustics, the fluid elements keep moving in the 
' direction set up by the initial condition. 

In order to proceed with a hydrodynamical description in which caustics do not form, a qualitative pressure gradient 
[l5| and thermal velocity scatter §,P^ in a collisionless matter have been discussed. Similarly, adhesion approximation 
(AA) [13 has been proposed, which is a model based on a nonlinear diffusion equation (i.e., Burgers's equation |18|). 
In AA, an artificial viscosity term is added to ZA; thus we can avoid caustics formation. The problem of structure 
formation has been discussed from the standpoint of AA 0, I20I I2H |2^ . where it is shown that the density 
divergence does not occur in AA and that a density distribution close to the N-body simulation can be produced. 
However, the origin of the viscosity in AA has not yet been clarified. 

Buchert and Dommguez |23( discussed the effect of velocity dispersion using the collisionless Boltzmann equa- 
tion [24| . They argued that models of a large-scale structure should be constructed by a flow describing the average 
motion of a multi-stream system. Then they showed that when the velocity dispersion is regarded as small and 
isotropic, it produces effective pressure or viscosity terms, and they consider that the isotropic velocity dispersion 
corresponds to the origin of the artificial viscosity in AA. Furthermore, in consideration of kinematic theory, they 
derived the relation between mass density p and pressure P, i.e., an equation of state. Buchert et al. [25| also showed 
how the viscosity term or the effective pressure of a fluid is generated, assuming that the peculiar acceleration is 
parallel to the peculiar velocity. Recently Buchert and Dommguez |26( provided an evaluation of the current status of 
adhesive gravitational clustering, which includes the above discussion, and they tried to improve past models. In their 
paper, they named their approach the Euler-Jeans-Newton (EJN) model. On the other hand, Domi'nguez p^ . |2S| 
proposed another approach. In these papers he clarified that a hydrodynamic formulation is obtained via a spatial 
coarse-graining in a many-body gravitating system, and that the viscosity term in AA can be derived by the expansion 
of coarse-grained equations. This model is named the Small-Size Expansion (SSE) model |26j. 

So far, with respect to the correspondence of the viscosity term with the effective pressure, and with regard to the 
extension of the Lagrangian description to various matter, the EJN model has been considered. Actually, Adler and 
Buchert [2^| have formulated the Lagran gian perturbation theory for a barotropic fluid. Morita and Tatekawa [30| . 
Tatekawa et al. [U, and Tatekawa [32l l33l| solved the Lagrangian perturbation equations for a polytropic fluid. 
However, it is still a open problem whether AA could realize the behavior of the EJN model in the proper way. It 
is known that the viscosity in AA decelerates the peculiar velocity in a dense region and avoids the formation of the 
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caustic, while the effective pressure in the EJN model also decelerates the motion of the fluid. But the fluid in the 
EJN model not only would decelerate but also might bounce in a dense region due to the effective pressure, which 
is known as the Jeans instability for cosmological fluid. In other words, we consider that AA would not realize the 
behavior of the EJN model. 

In this paper, to compare AA and the EJN model, we especially analyze the peculiar velocity in a cylindrical 
collapse model. Although we already compared AA with the EJN model with respect to density fluctuation in a 
past paper |34j . because we used an explicit method for solving partial differential equations, the accuracy of the 
numerical calculation in the dense region was not good. Therefore we could not observe the Jeans instability but, 
rather, the numerical instability. In this paper, instead of the explicit method, we apply the iterated Crank-Nicholson 
method [3{|. This method resolves the difficulty of solving the resulting implicit algebraic equations in the original 
Crank-Nicholson method and preserves the good stability properties. 

As our analyses show, although the peculiar velocity in AA decelerates due to the viscosity, the velocity does not 
oscillate. On the other hand, as a preliminary expectation, the peculiar velocity in the EJN model decelerates and 
also oscillates. Although the tendency of the evolution of the peculiar velocity would depend on the viscous parameter 
or the Jeans length, the behavior of AA is obviously different from that of the EJN model. We notice that while both 
AA and the EJN model certainly describe the quasi-nonlinear evolution well, the detail of the evolution is different. 
When we take a large value for the Jeans length, the fluid bounces to the outside. On the other hand, when we take 
a small value for the Jeans length, although the collapse of the cylindrical matter decelerates, a caustic might finally 
form at the center. Therefore it is problematic that the viscous term in AA is explained as an effect similar to the 
effective pressure term in the EJN model. For an explanation of the viscous term in AA, we will have to analyze more 
Lagrangian models. 

This paper is organized as follows. In Sec.^J we present Lagrangian perturbative solutions in the Einstein-de Sitter 
(E-dS) universe. First, we show linear perturbative solutions for dust fluid in Sec. Ill Al and in Sec. Ill Bl we mention 
the problem of ZA and show the solution of AA. Then we explain the EJN model in Sec. Ill Cl In Sec. lIIII we compare 
the evolution of the peculiar velocity in AA with that in the EJN model. Finally in Sec. IIVI we discuss our results 
and state our conclusions. 



II. THE LAGRANGIAN DESCRIPTION FOR THE COSMOLOGICAL FLUID 



In this section, we briefly present perturbative solutions in the Lagrangian description. In Newtonian cosmology, to 
introduce cosmic expansion, we adopt the coordinate transformation from physical coordinates r to comoving coor- 
dinates x, such as x = r/a(t), where a(t) is a scale factor. In Lagrangian hydrodynamics, the comoving coordinates 
x of the fluid elements are represented in terms of Lagrangian coordinates q as 

r = ax = a (q + s(q, t)) , (1) 

where s denotes the Lagrangian displacement vector due to the presence of inhomogeneities. With the Jacobian of 
the coordinate transformation from x to q, J = det(dxifdqj) — det(<5y + dsi/dqj), the mass density is described 
exactly as 

P = Pb J' 1 , (2) 



where ph means background average density. Furthermore we can decompose s into the longitudinal and the transverse 
modes, i.e., s = W q S + S T with V q • S T = 0. In this paper, we consider only the longitudinal mode for simplicity. 
The evolution equation for the longitudinal mode is written as follows [23 . l30| : 

V, • ( s + 2-s - j-7 Vx j ] = ^nGpbiJ- 1 - 1) , (3) 

where the dot above the variables denotes the partial derivative with respect to t. In general, it is very difficult to solve 
this equation for such reasons as the coordinate transformation or non-locality. In order to avoid this difficulty, we 
apply the perturbative approach and impose symmetry in Eq. ©■ Particularly in this paper we consider cylindrical 
symmetry, and the evolution equation in a cylindrical symmetric model is given in Appendix 1X1 
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A. The Lagrangian perturbation for dust fluid 

Zel'dovich derived a first-order solution of the longitudinal mode for dust fluid In the Friedmann Universe 
model, the solutions are formally written as follows: 

S(q, t) = D + {t)S+{q) + D_(t)S-(q) , (4) 

where D + (t) and -D_(i) mean the growing factor and the decaying factor, respectively. This first-order approximation 
is called the Zel'dovich approximation (ZA). In the case of the Einstein-de Sitter Universe, D+(i) and D_(t) are 
described as 

D+(t) oc e'\ (5) 
D_(t) cx r l . (6) 

Especially when we consider the plane-symmetric case, ZA gives exact solutions 0, . 



B. Adhesion approximation (AA) 

Although the Lagrangian approximation gives an accurate description until a quasi-linear regime develops, it cannot 
describe the model after the formation of caustics. After that, the nonlinear structure diffuses because the fluid 
elements keep moving in the direction set up by the initial condition. In order to avoid caustics formation, the 
adhesion approximation (AA) ^3] was proposed, which is a model based on a nonlinear diffusion equation (Burgers's 
equation 18]). In AA, an artificial viscosity term is added to ZA. In ZA, the equation for "peculiar velocity" u is 
written as follows: 



du 



+ (u-V x )u = 0, 



u = 



dx 
dDZ 



x 



(7) 

(8) 



where D + is the growing factor in ZA. To avoid caustics formation, in AA we add an artificial viscosity term to the 
right side of the Eq. Q, i.e., 



du 



(u ■ V x )u = i/Wiu . 



Now we introduce the Hopf-Cole transformation [sSEll such as 

u = -vV x (log (£,£>+)) ; 
then the Eq. © is changed to a diffusion equation: 

00 



Meanwhile, by using the inverse Hopf-Cole transformation, we obtain the solution of Eq. J5J: 



£>4 



D„ 



u(x{q,D'),D')dD' , 



Id 



3„/ fa-gQ 



G(x,x') 



fd 3 x'G(x,x>) 



G(x,x') — exp 



1 

2^ 



*o(x') + 



(x-x') 2 
2D 



(9) 
(10) 

(11) 

(12) 

(13) 
(14) 



where 



V x ^ {x) = s . 



(15) 
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We consider the case when the viscosity coefficient is quite small (y — > +0 (v ^ 0)). Within the limits of a small v, 
the analytic solution of Eq. © is given by 



u(x, t) 



y 



Eaia exp(-f^) 
where q a is the Lagrangian points that minimize the action 

I a = I(x,a;q a ) = S {q a ) + — - ^ = mm. 



2a 



det ( 8ij + 



d 2 S 
dqidq-j 



1-1/2 



S Q = S(q,t ), 



(16) 

(17) 

(18) 
(19) 



considered as a function of q for fixed x |2J, where the Roman character index denotes the Cartesian coordinate. 
In this paper, we consider the cylindrical symmetric case. For this case, we have to change the evolution equation, 
Eq. ©, slightly (see Sec. ITTTTjl . 



The EJN model 



Although AA seems a good model for avoiding the formation of caustics, the origin of the modification (or artificial 
viscosity) is not clarified. Buchert and Dommguez |23| argued that the effect of velocity dispersion is important in 
hindering caustics formation. They showed that when the velocity dispersion is still small and can be regarded as 
isotropic, it behaves as effective pressure or viscosity terms. Under the consideration of fluid kinematics, they proposed 
the effective equation of state as p oc p 5 ^ 3 . Also Buchert et al. [5^ showed how the viscosity term is generated by the 
effective pressure of a fluid under the assumption that the peculiar acceleration is parallel to the peculiar velocity. 
Moreover Buchert and Dommguez [2(| recently provided an evaluation of the current status of adhesive gravitational 
clustering, which is included in the above discussion. In their paper, they named their approach the Euler-Jeans- 
Newton (EJN) model. 

When we consider the polytropic equation of state P = up 1 , the first-order solutions for the longitudinal mode can 
be written in Fourier space as follows 30] ; for 7 ^ 4/3, 



S(K, a) cx a" 1 / 4 J ± 



I2C 2 



\K\ 



5/(8-67) 



C x |4-3 7 | 



(4-3 7 )/2 



(20) 



where S is the Fourier transformation of S, K is Lagrangian wavenumber, and J v denotes the Bessel function of 
order v, and for 7 = 4/3, 



S(K, a) oc a -i/ttvW^C4TO 



(21) 



where C\ = 47rG j Ob(ain) a in/3j C2 = K 7 /?b(ain) 7 a i n j ano - °i« means the scale factor given as an initial condition. 
When we take the limit k — > 0, these solutions agree with Eq. J3J. 

In this model, the behavior of the solutions strongly depends on the relation between the scale of fluctuation and 
the Jeans scale. Here we define the Jeans wavenumber as 



K, = 



/ 4irG Ph a 2 \ 
{dP/dp(p h )J 



1/2 



(22) 



The Jeans wavenumber, which gives a criterion for whether a density perturbation with a wavenumber will grow or 
decay with oscillation, depends on time in general. If the polytropic index 7 is smaller than 4/3, all modes become 
decaying modes and the fluctuation will disappear. On the other hand, if 7 > 4/3, all density perturbations will grow 
to collapse. In the case where 7 = 4/3, the growing and decaying modes coexist at all times. We can rewrite the 
first-order solution Eq. (|20|l with the Jeans wavenumber, i.e., 



S(K,a) cx a- 1 / 4 J. 



±5/(8-67) 



76 
|4 - 3 7 | K 3 



(23) 



5 



5 



-0 .005 




FIG. 1: Mexican-hat type model. The average of density fluctuation over the whole space becomes zero. 



In this paper, we analyze the first-order perturbation. When we consider cylindrical symmetric models, even if we 
deal with only the first-order perturbation, they cannot be analyzed in Fourier space, and we need to solve partial 
differential equations in real space with a numerical method (also see Appendix |A"|) . 



III. TIME EVOLUTION IN CYLINDRICAL MODEL 

For the cylindrical-symmetric case, dust collapse has been analyzed [l4|. Here we consider the evolution with ZA, 
AA, and the EJN models in the E-dS Universe. Hereafter, we define R and r as Eulerian and Lagrangian radial 
coordinates, respectively. At the initial time, we can identify the Lagrangian coordinate with the Eulerian one, i.e., 
R = r. Previous studies considered the collapse and/or evolution with the top-hat density distribution as the initial 
condition. Although the evolution of this model is easy to compute, the boundary condition becomes discontinuous. 
To avoid a discontinuity of the pressure gradient, we adopt the Mexican-hat type model (Fig. ^): 

8{R) =e(2- R 2 )e- R2/2 , (24) 

where R is Eulerian comoving radial coordinate. This model has several merits, for example the fluctuation is derived 
by the two times differential calculus of Gaussian, i.e., 

= e(2 - R 2 )e' R ^ 2 , (25) 
and the average of density fluctuation over the whole space becomes zero: 

I 2irR6(R)dR = . (26) 
J o 

For this model, from Eq. Q the solution of ZA is given as follows: 

S(a,r) = -aee- r2 ' 2 . (27) 

In our analysis, we set e = 1 due to an advantage of linear analysis, and the initial scale factor is set ao = 0.0167(= 
1/60) as the initial condition, where the initial density fluctuation at r = becomes S — 1/30 and the caustic appears 
at a — 1. On the other hand, for AA and the EJN model we can determine the initial longitudinal mode S and the 
initial "peculiar velocity" by preference as the initial conditions. In this paper they are made equal with those of the 
growing mode in ZA. Thus the initial conditions for AA and EJN model are given by 



S(a ,r) = -a ee- r ' 2 , (28) 



dS(a,r) 



da 



ee~ r I 2 . (29) 
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A. The adhesion approximation 



First we consider the evolution in the AA model. The evolution of the fluctuation is described by Eq. ©. For the 
cylindrical case, we slightly change the evolution equation. When we introduce cylindrical coordinates and assume 
the cylindrical symmetry, Burgers's equation is described as follows |39j : 



ud x u = 



~ 9 x (xd x u) - -4- 
X X 2 . 



1 



= vd x y~ d x(xu) 

where r and x are the time variable and the radial coordinate, respectively. Under a transformation such as 

2v 09 

u ~~Tdx~' 

Eq. H30|) is rewritten as 

d T = ~d x ( X d x 9) . 
X 



The generic solution for Eq. (|32|) is described by integral form: 



1 

2^ 



0(x,t) = fc(r)exp 
k{r) = 0(0, r), 
0(0, t) - 9 (x) = fc(0)ex P 



u(u, T)d,U! 



1 

2^ 



u(u>, r)doj 



(30) 

(31) 
(32) 

(33) 
(34) 
(35) 



In this paper we apply the cylindrical Burgers's equation (|3U|I for the adhesion model. In the E-dS Universe model, 
the evolution equation is given as 



du 
da 



du 



d 



U dR = V dR. 



R0R y ' 



Now we introduce Lagrangian displacement for the radial coordinate from Eq. Q : 

R = r + d r S{r, t) , 



(36) 



(37) 



where R and r mean Eulerian and Lagrangian radial coordinate, respectively, and S means Lagrangian displacement 
potential. Because we assume cylindrical symmetry and irrotational motion, the Lagrangian perturbation includes 
only the longitudinal mode. Also the "peculiar velocity" (Eq. (JSJ) is rewritten as 



u = 



0_ 

da 



dS_ 

Or 



(38) 



Using this "peculiar velocity," the evolution equation of the adhesion model is described by Eq. (|36() , and we examine 
the evolution of the "peculiar velocity." For comparison, we calculate cylindrical models with the viscous parameter 
v. Figure |21 shows the evolution of the "peculiar velocity" with v = 1/5, where we show the relation between Eulerian 
coordinates R and the "peculiar velocity" U. If we take a limit v — > 0, i.e, for the case of dust (ZA), the distribution 
of the "peculiar velocity" at any given time is the same as the initial distribution. On the other hand, from Fig.|2wc 
can see that the "peculiar velocity" decelerates monotonically due to the viscosity. 

We notice that Eq. I|36|l is described with Eulerian comoving coordinates. Thus it is complicated in a degree if we 
obtain the Lagrangian displacement (cf. Eq. In fact, for evolution of the fluctuation we need to consider the 

correspondence between Eulerian coordinates x and Lagrangian coordinates q for every grid at every time-step. 



B. The EJN model 



Next, we analyze evolution in the EJN model. If we consider a cylindrical model, because of mode-coupling in 
Laplacian, we cannot separate the perturbation to the spacial-dependent and the time-dependent term. Therefore we 
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FIG. 2: The evolution of the "peculiar velocity" in the AA model. Here we show the velocity distribution at a — 1/60 (initial, 
solid line) and a — 0.2, 0.4, • • • , 1.0. During evolution, the velocity decelerates because the viscosity resists the motion of the 
fluid. In AA, although the velocity decelerates, it does not oscillate. 



adopt numerical calculation for the evolution of the EJN model. In a previous paper |34j , we apply an explicit method 
for solving partial differential equations |4J|. However this method tends to produce numerical instability. To avoid 
the numerical instability, we adopt another method. Teukolsky |3^| discussed the iterated Crank-Nicholson method, 
which is one implicit scheme for numerical calculation. The method resolves the difficulty of solving the resulting 
implicit algebraic equations in the original Crank-Nicholson method and preserves the good stability properties. 
The "peculiar velocity" in comoving coordinates is described as 

v{a,q) = -^(V q S(a,q)) . (39) 
Because we consider a cylindrical symmetric model, we compute only radial velocity i.e., 

V r (a,r) = — y~L , (40) 



da \ dr 

where r is the Lagrangian radial coordinate. From the definition of the "peculiar velocity," we notice it does not change 
in the dust model (ZA). Because the solution in the dust model can be decomposed in time and spacial components, 
the "peculiar velocity" at any time is the same as the initial "peculiar velocity." Furthermore, the growing factor of 
the dust model is given by Eq. J5j. Therefore the longitudinal mode S in ZA monotonously increases (Eq. J57J)). 

Correspondingly we analyze the "peculiar velocity" in the EJN model. For the numerical calculation, we set the 
boundary condition at r = and r = 10 such as 



8S(a,r) 

S(a,W) 



dr 
dS(a,r) 



dr 



= 0, (41) 

r=0 

= 0. (42) 

r=10 



Because the distribution is asymptotically homogeneous, we set that the fluid does not move at the outside. The 
behavior of the EJN model strongly depends on the parameters k and 7. If k is very small, the solution is similar to 
that in ZA. On the other hand, if k is very large, the fluctuation oscillates and disappears. Following our previous 
method j3^|, we choose a reasonable value of the parameter 7, i.e., 7 = 4/3 and 5/3. In the case of 7 = 4/3, as we 
showed in Eq. 1|21[) . the solutions can be described simply and have both growing and decaying modes. In the case of 
7 = 5/3, Buchert and Domi'nguez j23| claimed that the isotropic velocity dispersion corresponds to the origin of the 
artificial viscosity in AA. Further, instead of k, we set the value of the Jeans wavenumbcr Kj (Eq. i|22|) ). 

Figures |3] and 0] show the "peculiar velocity" in the EJN model with linear approximation. We can see that the 
solution has both growing and decaying modes in the early stage, and that the position where the "peculiar velocity" 
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FIG. 3: The evolution of the "peculiar velocity" in the EJN model (7 = 4/3, Kj = 1.0). We show the time slices at 
a = 0.02, 0.04, 0.2. We also show the velocity distribution at a — 1/60 (initial, solid line). During time evolution, the 
peak of the "peculiar velocity" moves to the outside. Because of the effect of the pressure, the form of the distribution of the 
"peculiar velocity" changes from the initial one. The velocity decelerates and the falling fluid at the initial time bounces to the 
outside. 




r 



FIG. 4: The evolution of the "peculiar velocity" in the EJN model (7 = 5/3, Kj = 6.0). We show the time slices at 
a — 0.02, 0.04, 0.2. We also show the velocity distribution at a = 1/60 (initial, solid line). Similar to the case where 
7 = 4/3, during time evolution, the peak of the "peculiar velocity" moves to the outside. Because of the effect of the pressure, 
the distribution of the "peculiar velocity" oscillates at 0.02 < a < 0.2. 



is the fastest moves outward. Then the distribution of the "peculiar velocity" changes by the effect of the pressure. 
Moreover it is known that in Cartesian coordinates the solution for 7 > 4/3 always has a growing mode, and that 
the "peculiar velocity" is not zero at any time (especially in the case where 7 = 5/3 the perturbative solutions 
asymptotically approach those of the dust model). Therefore as seen Figs. and 0J the distribution of the "peculiar 
velocity" slightly oscillates, and the peak of the "peculiar velocity" moves to the outside at the early stage. Then, 
finally, caustics would form in the EJN model as in the case of the dust model (ZA). 
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FIG. 5: The evolution of the "peculiar velocity" in the EJN model (7 = 4/3, Kj = 5.0). We show the time slices at 
a — 0.2,0.4, • ■ ■ , 1.0. We also show the velocity distribution at a — 1/60 (initial, solid line). During time evolution, the peak 
of the "peculiar velocity" slightly moves to the outside. In this case, the "peculiar velocity" does not oscillate. However, the 
cylindrical dust would collapse, because the deceleration of the "peculiar velocity" is gentle. 

IV. DISCUSSION AND CONCLUDING REMARKS 

With respect to the distribution of the "peculiar velocity," we examine the correspondence between AA and the 
EJN model with cylindrical symmetry. In this analysis, even if we consider linear perturbation, it is hard to describe 
the solution with explicit form. Therefore we carried out numerical calculation, where to avoid numerical instability, 
we adopted the iterated Crank-Nicholson method. From our calculation, when we take small value for the Jeans 
wavenumber, the "peculiar velocity" in the EJN model oscillates due to the pressure. For the cosmological fluid, such 
oscillation is known as the Jeans instability. On the other hand, the "peculiar velocity" in AA is decelerated due to 
the viscosity, where there is no oscillation. Thus we can see that the behavior of AA is different from that of the EJN 
model. 

Furthermore, when we take a large value for the Jeans wavenumber in the EJN model, the evolution of the "peculiar 
velocity" is similar to that in AA (see Figure [SJ). However, it is well known that in Cartesian coordinates the linear 
perturbative solutions for 7 > 4/3 in the EJN model have a growing solution. Therefore the perturbation in the 
EJN model would eventually diverge. In other words, we can predict that even if we take a large value for the Jeans 
wavenumber in the EJN model, the radial motion of the fluid does not stop and finally collapses. Thus there also 
exists an essential difference between A A and the EJN model on this point, because it is known that the formation 
of caustic does not occur in AA. Hence AA cannot express the EJN model with propriety, and in order to clarify the 
origin of the viscosity term in AA, we should consider other effects besides isotropic velocity dispersion or isotropic 
effective pressure. 

Recently, Buchert and Dommguez p6| discussed adhesive gravitational clustering and tried to provide a clear 
explanation of the assumuptions for adhesion approximation. They applied the Eulerian and Lagrangian expansions 
to the non-perturbative regime and proposed a new non-perturbative approximation. When we analyze this new 
approach with both analytic and numerical methods, we may be able to explicate the origin of the artificial viscosity 
in AA. Also, as a future work, we would describe the evolution of the density fluctuation in a nonlinear regime with 
a semi-analytic method. 
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In Sec. Ill CI we noticed that Eq. © is hard to solve because of the coordinate transformation or non-locality. In 
this appendix we rewrite the equation with cylindrical symmetry. We can introduce Lagrangian displacement for the 
radial coordinate by Eq. (|37|) . The spacial derivative is rewritten with the Lagrangian coordinate, such as 



APPENDIX A: THE EVOLUTION EQUATION IN CYLINDRICAL MODEL 



d_ 

Or 



8Rd_ 
~dr~ ~dR 



(1 + d 2 r S) 



d_ 
dR' 



(Al) 



Therefore, the derivative is changed as 



d_ _ 1 d 
dR ~ l + d^Sdr' 



(A2) 
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The divergence of the "peculiar velocity" with a Eulerian coordinate becomes a little complicated: 

1 d 



V * U = RdR {RdRU) 
1 1 



d 



{{r + d r S)d r u} 



(r + d r S) 1 + d?S dr 

We decompose the Jacobian of the coordinate transformation to the order of the perturbation. 

,7=1 + + J(2) + jO) 



where 



J (1) = d?S = d?S+-d r S, 

T 

J (2) = \ Ms) (ap) - (d^s) (djdiS)] 



= - (d r s) (d 2 r s) , 



J 



(3) _ 



det (didjS) = 0. 



(A3) 

(A4) 
(A5) 

(A6) 
(A7) 



Most of above changes affect only higher-order approximation. In the EJN model, we must consider the change of 
Laplacian in Lagrangian space, i.e., 



r dr 
1 
/• 



df(r) 



dr 

(f'(r)+rf"(r)) , 



(A8) 



where the prime denotes the partial derivative with respect to r. Using the above deformations, we obtain the evolution 
equation for the EJN model with cylindrical symmetry. The linear perturbative equation for the EJN model is written 
as 



ld_ 

r dr 



r{S' + 2-5' 

a 



j2 dr \r dr 



= -4»G ft -l|-[rS'] . 
r or 



When we choose an appropriate boundary condition, we can rewrite Eq. I|A9|) as 



a a z 



5" + -S' 

r 



-4irGp b S . 



(A9) 



(A10) 



When we treat the EJN model in Cartesian coordinates, we can apply a Fourier transformation to the perturbative 
equations easily. However, with the cylindrical symmetry, because the fourth term on the left-hand side in Eq. (|A10|) 
derives convolution, even if we take only linear perturbation, we cannot obtain analytic solutions and we have to solve 
by numerical calculation. 



